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The large mass of the ninth pseudoscalar meson, the rj , is believed to arise from the combined 
effects of the axial anomaly and the gauge field topology present in QCD. We report a realistic, 
2+1 flavor, lattice QCD calculation of the rj and r/' masses and mixing which confirms this picture. 
The physical eigenstates show small octet-singlet mixing with a mixing angle of 6 — —14.1(2.8)°. 
Extrapolation to physical light quark mass gives, with statistical errors only, m, = 573(6) MeV and 
m n > = 947(142) MeV, consistent with the experimental values of 548 MeV and 958 MeV. 

PACS numbers: 11.15.Ha, 11.30.Rd, 12.38.Gc 14.40.Be 



The relatively large mass of the ninth pseudo-scalar 
meson, the 77', provides a significant challenge for quan- 
tum chromodynamics (QCD), the component of the stan- 
dard model which describes the interactions of quarks 
and gluons. On a naive classical level there are nine con- 
served axial currents. Given the vacuum breaking of the 
symmetries which these currents generate, this should 
imply the existence of nine Goldstone bosons, a conclu- 
sion inconsistent with the large splitting between the 8 
octet mesons, n^, tt°, K , K°, K , r\ and the singlet 
ff Q- Unique among these nine currents, the U(l) axial 
current, corresponding to the singlet 77' meson, has an 
anomalous divergence at the quantum level. However, 
to arbitrary order in perturbation theory this anoma- 
lous divergence vanishes at zero momentum, continuing 
to imply that the masses of all nine pseudoscalar mesons 
should vanish in the limit of vanishing quark mass. It 
is only with the discovery of instanton configurations 
with non-trivial topology Q that a mechanism Q be- 
came available that could explain the large 77' mass. 

While these important developments suggest possible 
consistency between QCD and the value of the 77' mass 
a direct demonstration of the required anomaly-driven, 
octet-singlet splitting has been lacking. In this paper we 
present the first such demonstration in the realistic case 
of three light dynamical quarks. 

The critical role of disconnected diagrams in the study 
of the 77 and 77' and the severe difficulties they introduce 
have been recognized for more than 15 years [1,0]. Posi- 
tivity requires the quark propagators that appear in the 
connected diagrams to decrease exponentially with in- 
creasing time separation. For mesons this fall-off roughly 
matches the exponential time dependence of the massive, 
Euclidean-space meson propagator and good numerical 
signals can be seen over a large range of times. For terms 



in which the source and sink of the meson propagator are 
not joined by quark propagators, the needed exponential 
decrease comes from increasingly large statistical cance- 
lations implying a rapidly vanishing signal-to-noise ratio. 
These difficulties have impeded earlier work on this 
topic which has employed indirect methods or not exam- 
ined the physical case of up, down and strange dynamical 
quarks; see also Ref. [lOj . 



SIMULATION DETAILS 

Our calculation uses the Iwasaki gauge and domain 
wall fermion actions, a 16 3 x 32 space-time volume with a 
fifth-dimensional extent of 16 and /3 = 2.13. giving an in- 
verse lattice spacing 1/a = 1.73(3) GeV [11(. We analyze 
three ensembles of gauge configurations with light sea 
quark mass m t = 0.01, 0.02, 0.03 12]. (All dimensionful 
quantities are given in lattice units except when physical 
units are declared.) These values of m; yield pion masses 
of 421, 561 and 672 MeV, respectively. The 0.01 and 
0.02 ensembles were generated using the physical strange 
quark mass m s = 0.032 [13T ]. The mi = 0.03 ensemble 
was reported as RHMC II in Ref. [ll| with m s = 0.04. 
For this ensemble we use reweighting to change m s from 
0.04 to 0.032 in 20 mass steps [L^. 

We use a Coulomb gauge fixed wall source and sink 
for the quark propagators. Because of the difficulty of 
computing the disconnected graphs, large statistics are 
required. Therefore, we calculate propagators for sources 
on each of our 32 time slices. The large number of Dirac 
operator inversions (32 x 12) that must be performed on 
a single gauge configuration is accelerated by computing 
the Dirac eigenvectors with the smallest 35 (mi = 0.01) 
or 25 (mi = 0.02, 0.03) eigenvalues and limiting the 
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FIG. 1: Five diagrams appearing in the 77 and 77' correlation 
functions. They are Cu(t), C 3S {t), T> u {t), V S3 (t), and V ts (t) 
respectively from left to right and top to bottom. The solid 
lines are quark propagators and the solid circles 75 insertions. 



conjugate gradient inversion to the remaining orthogonal 
subspace. This results in a 60% speed-up for mi = 0.01. 
We study 300 configurations separated by 10 molecular 
dynamics time units for mi = 0.01 and 0.02, and 150 
configurations separated by 20 time units for mi = 0.03. 

We compute four Euclidean space correlation functions 
between two pseudoscalar operators Oi and O s : 

1 31 

C(t) a p = — ^2(O a (t + a,pe{l,s}, (1) 

i'=0 

summed over the 32 source locations. Here O s = sj^s 
and Oi = (U75U + d^d)/^, both SU(2) singlets. 

The matrix C(t) can be expressed in terms of the five 
amplitudes represented by the diagrams shown in Fig. [T] 
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This equation shows that neither Oi nor O s creates an 
energy eigenstate of QCD. They mix with each other 
through the disconnected diagram V s i = T>i s . The usual 
expectation that such disconnected graphs are small does 
not apply here. Figure [2] shows these amplitudes ver- 
sus time for the mi = 0.01 ensemble. The disconnected 
graphs decrease more slowly than the connected graphs, 
changing the pattern of SU(3) flavor symmetry breaking. 

Inserting a sum over states into Eq. [1] and assuming 
this sum is dominated by the 77 and rj' for large t we 
obtain 



C(t) = A T D{t)A. 
where the overlap matrix A is given by: 



A 



(n\Oi\o) ( v \o s \o) 

( V '\Oi\0) (rf\O,\0) 



(3) 



(4) 



and 



and D(t) is a diagonal matrix with elements e~ 
^"V*. We chose A real, possible because C{t) is real. 

Now define a second operator basis with definite SU(3) 
properties: the octet Os = (uj 5 u + dj 5 d — 2s7 5 s)/\/6 
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FIG. 2: Results for the five contractions which enter the 77-77' 
correlator calculated using the mi = 0.01 ensemble. 



1 



-e — © — © — © — $ — $ — ^ — o 



■q 1 — e- 

Tl' 1 — 

0.8 
0.6 
0.4 
0.2 


2 4 6 8 10 12 14 
t 

FIG. 3: Effective mass plot for the 77 and 77' states from the 
m; = 0.01 ensemble. 



and the singlet 0\ — (u^f 5 u + d^ 5 d + s7 5 s)/v / 3- We 
will use the Roman indices a and b, for these operators, 
e.g. {O a } a= s,i to distinguish them from the earlier basis 
{O a } a= i tS . Equations analogous to Eqs. [TJ [3] and H] will 
be obeyed if this second basis with a, fee {8,1} is used. 

We can determine the two masses and the four real 
elements of the matrix A by fitting our data to Eq. [3] 
over an appropriate range of time t. To determine this 
range we examine the product: 



C(i ) _1 C(i) = A~ 1 D(t - t ) A, 



(5) 



implying C(to)~ 1 C(t) is similar to a diagonal matrix 
whose eigenvalues are exponentials of the masses of in- 
terest. We find the best results if t — t is large, giving a 
clean separation of the larger, more accurate r\ eigenvalue 
and the smaller eigenvalue associated with the noisy 77'. 
Figure [3] shows the eigenvalues obtained from Eq.[5j Here 
we plot the logarithm of the ratio of each eigenvalue eval- 
uated at t and t + 1 with to = 2. The choice to = 2 and 
3 < t < 7 gives a recognizable plateau for m v and m. q > . 
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r? — 77' MIXING 

It is customary to treat the physical r\ and 77' states 
as mixtures of the pseudo-scalar octet and singlet states 
which appear in the SU(3) symmetric limit and to in- 
troduce an angle which specifies this mixing. In the 
present calculation we can examine the validity of this 
mixing model and attempt to determine 0. Consider the 
SU(3) symmetric limit mi — m s and let |8) sym and |l) S ym 
be these lowest energy octet and singlet states with en- 
ergies Eg and E\. We justify this mixing model by as- 
suming that when mi ^ m s the only important effects 
are a subset of those implied by first-order perturbation 
theory: first-order energy shifts and first-order mixing of 
states but only for those cases enhanced by the relatively 
small energy denominator E\ — Eg. To zeroth order in 

1/2 

m s — mi we can write S ym(a|Ofc|0) = Z a S a b and we as- 
sume this relation is unchanged by the first order effects 
of m s — mi on the vacuum state — again neglecting mix- 
ing not enhanced by the factor 1/(E\ — Eg). 
These assumptions imply that 



cos(0) 
sin(0) 



- sin(0) 
cos(0) 




(6) 



and that the overlap matrix A can be written: 



A = 



Zg cos(0) -Z{ sin(0) 
z\ sin(0) z\ cos(0) 



(7) 



for A in the Og - 0\ basis. The columns of A are thus 
orthogonal and, if Og and 0± are normalized by multipli- 
cation by Z s and Z 1 X ^ 2 , the resulting overlap matrix 
A will be orthogonal. Using the results below, we find 
for the dot product between the columns of A: -0.016(9) 
and -0.012(4) for the mi = 0.01 and 0.02 ensembles. 

We can also extract an effective mixing angle 0(i) from 
Eq. [5l This equation determines each row of A up to 
an arbitrary constant. However, these two undetermined 
normalization factors as well as the factors Z^ 2 and Z\^ 2 
cancel from the product A^A^g/ A^gA^i, a combination 
which equals — tan 2 (0). The resulting angle is shown in 
Fig. |U The small value of in the Og and 0\ basis 
demonstrates the large role played by the disconnected 
diagrams. Had we omitted the disconnected diagrams, 
the matrix A would have been diagonal in the Oi and 
O s basis giving sin(0) = —y/2/3 or = —54.7°, very 
different from our = -14.1(2.8)°. 



FITTING RESULTS 

We fit our four correlation functions C a b(t) in two 
steps. First, using 3 < t < 7 we determine the two masses 
and the four elements of A. Second, we fix 



15 
10 
5 


-10 
-15 
-20 
-25 



ii " 



2 4 6 

t 



10 12 



FIG. 4: The 77 — 77' mixing angle 9(t) determined from Eq.[5]for 
the mi = 0.01 ensemble. While the errors are large, the data 
is consistent with a single value of about -10° for 3 < t < 7. 
(Note, 8(t) is undefined at to = 2 and off scale at t = 0.) 



TABLE I: Meson masses for the mi/m s = 0.03/0.04 ensemble 
and at the reweighted value m B s ca — 0.032 for two values of 
the valence strange quark mass mj al = 0.03 and 0.04. Here 
and below only jackknife, statistical errors are given. 







mj al 


m n 


77V 


0.04 


0.3907(9) 


0.03 


0.3907(9) 


0.716(49) 




0.04 


0.4316(16) 


0.713(67) 


0.032 


0.3899(11) 


0.03 


0.3899(11) 


0.688(60) 




0.04 


0.4328(20) 


0.694(126) 



A to that determined in the first step and fit the 7777 ele- 
ment of the transformed matrix \{A T )~ 1 C(t)A~ 1 ] m over 
the larger range 5 < t < 15 to determine more accurately 
m n . For each fit we minimize \ 2 computed from the full 
covariance matrix, which includes the statistical correla- 
tions between each measured propagator at each of the 
time separations used. We treat each configuration as 
independent but check for autocorrelations by grouping 
the data into blocks of size up to 10 and find consistent 
errors. As a test for long autocorrelations, we compare 
the first and second halves of our data and find consistent 
results. Using random sources, Refs. (s , 15 1 suggest these 
disconnected correlators show large statistical excursions. 
We do not see this behavior. Our standard wall sources 
give disconnected and connected propagators which fol- 
low similar, properly-sampled distributions. 

For the mi — 0.03 ensemble, we reweight the correla- 
tion functions to change m| ea from 0.04 to 0.032 and list 
the results in Tab. fl] We then linearly interpolate the 



m v and m v 



resulting m^ and m^ with strange valence quark masses 
of 0.03 and 0.04 to the point mf = 0.032. Table |H] lists 
the resulting masses for the octet states, n, K, 77, and 
the singlet state 77' for each ensemble. The final column 
shows m v determined by the Gcll-Mann-Okubo(GMO) 
formula Zrrvz + m\ — Am 2 K using our values for and 
to if . The good agreement with this first order formula is 
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TABLE II: Masses in lattice units for the nonet of pseudoscalar mesons. 



m;(conf) 




m K 


m n 
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m„(GMO) 


0.01(300) 
0.02(300) 
0.03(150) 


0.2441(7) 
0.3251(6) 
0.3899(11) 


0.3272(7) 
0.3633(6) 


0.3572(24) 
0.3787(11) 
0.3988(13) 


0.600(45) 
0.605(36) 
0.689(73) 


-8.3(2.6)° 
-5.5(1.4°) 


0.3505(10) 
0.3752(9) 



consistent with our small octet-singlet mixing. 




FIG. 5: Extrapolation of m v , m n i (upper) and 8 (lower) to 
physical light quark mass (and a negative input mass mi). 

In Fig.[S]we show a linear extrapolation of m v > and 
as a function of to the physical value of , consistent 
with NLO chiral perturbation theory. (Note, the curva- 
ture of the m, ; fit is barely visible.) We find m v — 573(6) 
MeV and m^i — 947(142) MeV, where the errors are 
statistical. To verify our choice of m s , we extrapolate 
the kaon mass and find the physically consistent value 
497.4(7) MeV. Also shown is a similar linear extrapola- 
tion for 8 giving 8 — —14.1(2.8)°, in agreement with the 
range —10° to —20° of phenomenological values (l6| . 

We have described a 2+1 flavor calculation of the 
masses and mixing for the r\ and rj mesons finding re- 
sults agreeing within their 15% error with experiment. 
The near orthogonality of the mixing matrix A is con- 
sistent with physical states which are simple mixtures 
of SU(3) octet and singlet states. Given our large sta- 
tistical errors we have not analyzed the smaller system- 



atic errors arising from our single lattice spacing, large 
light quark masses and finite volume which other calcula- 
tions [Ill|l3,ll2| suggest are « 4%, 5% and 1%. However, 
to this accuracy our calculation demonstrates that QCD 
can explain the large mass of the ninth pseudoscalar me- 
son and its small mixing with the SU(3) octet state. 
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